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We propose a model of dark matter identified with a pseudo-Nambu-Goldstone boson in the 
dynamical supersymmetry breaking sector in a gauge mediation scenario. The dark matter particles 
annihilate via a below-threshold narrow resonance into a pair of R-axions each of which subsequently 
■ decays into a pair of light leptons. The Breit-Wigner enhancement explains the excess electron and 

positron fluxes reported in the recent cosmic ray experiments PAMELA, ATIC and PPB-BETS 
without postulating an overdensity in halo, and the limit on anti-proton flux from PAMELA is 
£NJ 1 naturally evaded. 
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I. INTRODUCTION 



Existence of non-baryonic dark matter as the dominant component of matter in the universe has been established 
by numerous observations. The origin of the dark matter, however, has not been identified yet, and its nature is 
arguably the most important problem in particle and astrophysics. 

Recent observations of ATIC 1] and PPB-BETS [2| balloon experiments show the existence of a bump in a 300- 
800 GeV energy region of e~ + e + flux in cosmic ray. The interesting astrophysical possibilities for the origin of 
Oh the excesses are nearby pulsars [H, [U, H| or supernovae remnants @ . The most exciting interpretation of the excess, 
however, is the annihilation and/or decay of the dark matter with a mass in a TeV range. It is remarkable that it 
explains simultaneously the anomalous excess of e + flux in PAMELA experiments Q • 

From the theoretical point of view, it is very interesting to explain the dark matter in the supersymmetric (SUSY) 
extension of the standard model (SSM), since the SUSY models naturally possess two types of candidates for the dark 
matter. One that has been discussed extensively in the literature is the lightest SUSY particle called as the LSP [8[, 
and the other often overlooked is stable composite "baryons" in a dynamical supersymmetry breaking sector [9( . The 
former case, however, implies that the masses of the gluino and squarks are much larger than a TeV range, which 
^vqj . causes serious problems in discoveries of SUSY particles at the LHC. The latter case predicts most likely the mass of 
the dark matter to be in at least about 30 TeV, and hence, it seems difficult to explain the dark matter with a TeV 
. mass. 

' In this paper, we propose a model where the dark matter with a TeV mass is nothing but pseudo-Nambu-Goldstone 
bosons generated by the SUSY breaking sector (as in the latter case above) . Surprisingly, this model provides solutions 
1 to the two puzzles in the recent cosmic ray experiments. The first puzzle is that the required annihilation cross section 
in the galactic halo is much larger (by a factor of 0(100)) than the one appropriate to explain the dark matter relic 
density precisely measured by the WMAP experiment . The second puzzle is that the PAMELA experiment sees 
5— i ' no excess in the anti-proton flux while it sees an excess of the anti-electron flux. 



As we will see, in the Nambu-Goldstone dark matter scenario, the observed dark matter abundance is achieved 
only if the annihilation process occurs near the pole of a narrow resonance. This inevitably evokes the Breit-Wigner 
enhancement of the dark matter annihilation [TjJ (see Refs. [IH, [l3j for earlier attempts) which explains the so-called 
boost factor. Furthermore, we will see that the dominant final state of the near-pole annihilation is a pair of the 
R-axions each of which subsequently decays into a light lepton pair. Therefore, this model also provides a concrete 
example of the scenario flU [Ta ] explaining the second puzzle in the PAMELA data. 

The organization of the paper is as follows. In scctionHTl we will discuss generic features of the dark matter in 
the SUSY breaking sector in the light of the recent cosmic ray experiments. In section lLUl we propose the Nambu- 
Goldstone dark matter where the dark matter annihilates via a narrow resonance into an R-axion pair. The final 
section is devoted to our conclusion. 



II. GENERAL DISCUSSION ON THE HIDDEN SECTOR DARK MATTER 

The idea of the dark matter in the SUSY breaking sector in gauge mediation models was first sketched in Ref. [|[. 
When the SUSY breaking sector possesses a global symmetry, the lightest particle which is charged under the global 
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symmetry is stable and can be a candidate of the dark matter. In earlier attempts, the mass of the dark matter was 
postulated to be of the order of the dynamical SUSY breaking scale, i.e., around 30TeV, which is the lowest possible 
scale realized in gauge mediation scenarios [IB, E3, 0, EH. However, such a heavy dark matter is not favorable to 
explain the observed bump in a 300-800 GeV energy region of e~ + e + flux in cosmic ray [j], Therefore, in order 
to obtain a viable dark matter model in the SUSY breaking sector, we need to consider some mechanisms to realize 
a "light" dark matter candidates in the hidden sector. 

The simplest possibility to obtain such a light particle is to introduce a small coupling so that the dark matter 
candidates have a small mass, i.e., 

todm = eAsusy, (1) 

where m^M is the mass of the dark matter which is suppressed by a small coupling e compared with the SUSY 
breaking scale Asusy- The dark matter with mass in a TeV range requires 

^lo-a- 2 ), (2) 

when the SUSY breaking scale is around 30 TeV. However, the introduction of the small coupling naively ends up with 
a too small annihilation cross section to explain the observed dark matter abundance. That is, the naive estimation 
of the annihilation cross section of the dark matter at the freeze-out time, 
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is too small to explain the observed dark matter abundance which requires, 

av Icl ~ W- 9 GeV~ 2 . (4) 

Here, we are assuming that the final state is lighter particles in the hidden sector which eventually decay into the 
SSM particles. Besides, since we are assuming the models with gauge mediation, the coupling between the hidden 
sector and the SSM sector is rather suppressed (see also discussion in section HTT|) . 

A more ambitious possibility is to identify the light dark matter to the pseudo-Nambu-Goldstone bosons resulting 
from a spontaneous breaking of an approximate global symmetry in the SUSY breaking sector in analogy with the 
pions in QCD. In this case, we do not need to introduce small couplings to realize the light dark matter. However, 
the naive estimation of the annihilation cross section of the dark matter is again suppressed, i.e., 

1 / todm V 1 

O-Vrcl ~ — T —2—, 5 

16?r V a susy/ mf, M 

where we have assumed the breaking scale of the approximate global symmetry to be of the order of Asusy- 

Theses lessons tell us that the annihilation cross section of the dark matter must be enhanced than the above 
naive expectations. As an interesting possibility, such enhancement can be realized if we assume that the dark 
matter annihilates via a narrow resonance with mass M ~ 2?tidm- This observation, in turn, suggests the possible 
enhancement of the dark matter cross section in the galactic halo by the Breit-Wigner enhancement mechanism [ll[ . 
In the following section, we construct a model of the Nambu-Goldstone dark matter where these possibilities are 
realized. 



III. NAMBU-GOLDSTONE DARK MATTER 



In this section, we construct an explicit model of the Nambu-Goldstone dark matter based on a dynamical SUSY 
breaking model. For that purpose, we consider a vector-like SUSY breaking model developed in Ref. [2(|. As we will 
see, this model possesses all the necessary ingredients to realize the Nambu-Goldstone dark matter model where the 
Breit-Wigner enhancement explains the effective boost factor and the R-axion final state explains the no excess in 
anti-proton flux. 
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A. Vector-like SUSY Breaking Model 

The vector-like SUSY breaking model is based on an SU(2) gauge theory with four fundamental representation 
fields Qi(i = 1, • • • ,4) and six singlet fields Sij = —Sji = 1, ■ ■ • ,4) [20]. In this model, the SUSY is dynamically 
broken when the Q's and 5's couple in the superpotential, 

W = XijSijQiQj, (i < j), (6) 

where Ay denotes coupling constants. The maximal global symmetry this model may have is 5P(4) ~ 50(6) symmetry 
which requires Ay = A. The SUSY is broken as a result of the tension between the F-term conditions of 5's and Q's. 
That is, the F-tevm conditions of Sij, dW/dSij — XijQiQj — 0, contradict with the quantum modified constraint 
Pf(My-) = A d where My- denote composite gauge singlets made from QiQj. Especially, when the coupling constants 
Xij are smaller than unity, the SUSY is mainly broken by the F-term. of a linear combination of the singlets Sy . 

The effective theory below the dynamical scale Ad yn is well-described by the gauge singlets My and Sij with the 
effective superpotential, 

W cS = XijA dyn SijMij + X (Pf(M) - A 2 dyn ) , (i < j), 

= *AAdynS A M A + X I ^l-Adyn), (7) 

A=0-5 VA=0— 5 / 

where A" is a Lagrange multiplier field which enforces the quantum modified constraint, and we have rearranged the 
Sij and Mij by using appropriate linear combinations in the last expression. Here, we have assumed that the effective 
composite operators Ma are canonically normalized (up to order one ambiguity in the coefficient that we will neglect 
in the following). 1 

Now let us assume that the SUSY breaking sector possesses an 50(5) C 50(6) global symmetry, and take A = Ao 
and A' = A a= i_5 with A < A'. In this case, the lightest particle which is charged under the 50(5) symmetry is stable 
and can be the dark matter candidate. Under these assumptions, the quantum modified constraint is solved by, 



Mo = jA 2 dyn - £ Ml (8) 



y 0=1—5 

By plugging Mq into the effective superpotential in Eq. (|7|). we obtain 



W cS ~ XA 2 dyn S - t; s oMZ+ £ A' A dyn S a M a + 0(M%). (9) 

a— 1 — 5 a— 1 — 5 

Thus, in terms of the low energy effective theory, the SUSY breaking vacuum is given by, 

F So = XA 2 dyn , S a = 0, M a = 0. (10) 

B. Dark Matter without R-symmetry Breaking 

1. R- symmetric spectrum 0} the light particles 

Before introducing the R-symmetry breaking, it is worth considering the model with no R-symmetry breaking, i.e., 
(Sq) — 0, which clarifies the necessity of the R-symmetry breaking. In the case of the R-symmetric vacuum, the mass 
spectrum is given as follows. First, the lightest particle which is charged under the 50(5) comes from the scalar 
components of M a whose masses squared are given by, 



in 



,2 



± = (A' 2 ±A 2 )A^ yn , (11) 



where the minus sign corresponds to the real component of the M a scalar. The A dependence comes from the SUSY 
breaking effect coming through the SqM 2 coupling in Eq. ([9]). On the other hand, the scalar part S a does not receive 
the SUSY breaking effects, and has the same mass with the fermion components of 5 Q and M a , i.e., 

m Sa = m Ma = A' A dyn . (12) 



1 If we use the naive dimensional analysis 12 ll . Ad yn is replaced with A^ yn /An without affecting the following discussions. 
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Therefore, we find that the dark matter is given by Re[Mj. 

The masses of the Sq components require attention. Since the scalar component corresponds to a classical flat 
direction, its mass vanishes at the tree-level. The one-loop Coleman- Weinberg potential of Sq, however, gives rise to 
the mass of So as (see appendix^]) , 

A 3 

ms ° ~ J^jx Adyn - (13) 

In contrast, the fermion component of Sq contains the goldstino which acquires a very small gravitino mass by coupling 
to supergravity. 

By putting it all together, we find that the masses of the dark matter as well as the other components of S a , M a 
and Sq are parametrically lighter than the dynamical scale A<j yn . For example, we obtain a light dark matter for small 
couplings, 

m DM ^eA SUSY , e = 0(A 1/2 , A' 1/2 ), (14) 

where we have used AgusY = A 1/,2 Ad yn and A < A'. Thus, the dark matter with a mass in a TeV range can be achieved 
for 

e = 10-( 1 - 2 ), (A,A' = 10-( 2 - 4 >). (15) 

The above spectrum also poses the other possibility discussed in the previous section, i.e., the pseudo-Nambu 
Goldstone boson dark matter. We can see it by taking the limit of 50(6) global symmetry, i.e., A — > A'. There, the 
mass of the dark matter in Eq. ([TTj) vanishes. This shows that the dark matter is nothing but the pseudo-Nambu- 
Goldstone boson of the spontaneous breaking of 50(6) app — > SO(5) with a breaking scale Ad yn (see Eq. (JHJ). In this 
case, we obtain the dark matter with a TeV mass for 

A'- A = 0(l(T (2 - 4) ), (16) 

while keeping A and A' of the order of one. 



2. Dark matter annihilation without resonance 



Now, let us consider the dark matter annihilation. For, nis < m-DM, the dark matter Re[M a ] dominantly annihilates 
into So scalar via the F-term potential \ms a S a — A5oM a | 2 (see Eq. ©). The amplitude of this process is given by, 

M = A 2 + A 2 ms % = A 2 ^-s-, (17) 

t — m s t — nij 

where t denotes the momentum transfer. The first term comes from the four-point interaction and the second term 
from the t-channel exchange of the S a scalars. In the 5-wave limit, the momentum transfer is given by, 



t = -ml M (3j, l3 f = Jl-^s.~l, (18) 

V m DM 



and the cross section is given by, 



/3f v rei 

0-«reI = 7^-^77; WTT* 



87r2(2m D M) 2 A \t - m 2 Sa 

f I m DM ^ 



a 4 m 



32tt m% M 



where the final approximation is valid for todm ~ ms, 



(19) 



2 For mg Q > mrjMi the dark matter dominantly annihilates into the gravitinos with a much more suppressed annihilation cross section. 
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Prom this expression, we confirm that the cross section of the "light" dark matter, i.e. A, A' = 10~( 2 ~ 4 ) (Eq. ([15)) ') 
or m DM ~ 1 TeV and mg a ~ 30TeV (Eq. p^| ). is highly suppressed. So, we need to look for an appropriate narrow 
resonance so that the cross section is sufficiently enhanced. Interestingly, in the case of the Nambu-Goldstone dark 
matter, there is a candidate for such a resonance, the scalar part of So- The Eq. (fT3"|) shows that the mass of Sq scalar 
can be also in a TeV range for A < 1, and hence, the So mass can satisfy ms — 2 todm with a careful tuning. Thus, if 
the dark matter annihilates via the So resonance, the annihilation cross section can be drastically enhanced from the 
one given above. However, for this process, the R-symmetry must be broken, since the R-charge of So is 2, while that 
of M a is 0. Motivated by these observations, we will extend our analysis to the model with the R-symmetry breaking. 



C. Nambu— Goldstone Dark Matter with R-symmetry breaking 

1. R-symmetry breaking 

Now, let us consider spontaneous R-symmetry breaking. For simplicity, we assume that the R-symmetry is broken 
by effects of higher dimensional operators of So in the Kahler potential, 

* = ' S °' 2 + W-W + '"' (20) 

where A's denote the dimensionful parameters and the ellipsis denotes the higher dimensional terms of So- The 
positivity of the coefficient of the quartic term is crucial to destabilize the R-symmetric vacuum at So = 0. Notice 
that the above Kahler potential provides an effective description of a quite general class of the models with spontaneous 
breaking of the R-symmetry breaking. Especially, when the above Kahler potential results from radiative corrections 
from physics at the scale Ad yn , the dimensionful parameters are expected to be, 



J_ A L_ J_ A L_ rm 

A! le^A^' a| le^AV 1 > 

where dimensionless coefficients ca.q are of the order of unity. In appendixlBl we demonstrate an explicit perturbative 
model which breaks the R-symmetry in a similar way studied in Ref. [22j ■ 

From the above Kahler potential, the R-symmetry is spontaneously broken by the vacuum expectation value of the 
scalar component of So; 

< s °> - = ^>>-< = >• < 22 ' 

where we have introduced the R-symmetry breaking scale fn — 0(Ad yn ) and define the R-symmetry so that (So) > 0. 
At this vacuum, the scalar component of So is decomposed into a flaton s and the R-axion a by, 

S = ^=(f R + s)e ia ^. (23) 



Then, the mass of the flaton is given by, 



AA?„„A? _AAjL 



m. 



where we have used F$ = AA^ yn and assumed Eq. (|2"Tj) with a — c 6 = 0(1). Therefore, the flaton can be in a TeV 
range for A ~ 1 and C4 ~ 1, which is a crucial property for the flaton to make the narrow resonance appropriate for 
the dark matter annihilation. 

On the other hand, the R-axion mass is much more suppressed and mainly comes from the constant term in the 
superpotential which breaks the R-symmetry explicitly. 3 In the supergravity with (almost) vanishing cosmological 



In this study, we assume that the messenger sector of the gauge mediation also respects the R-symmetry. Otherwise, the radiative 
correction to the Kahler potential of So from the messenger sector gives rise to the dominant contribution to the R-axion mass. The 
R-breaking mass from the Higgs sector, on the other hand, is smaller than the one in Eq. (1251) . even if the so-called /i-term does not 
respect the R-symmetry. 
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constant, the R-axion acquires a small mass [23 



ml^^f^- (25) 

JR 

In the case of the low-scale gauge mediation with the dynamical SUSY breaking scale around 30 TcV, the R-axion 
mass is tens to hundreds of MeV range. 



2. Spectrum and interactions of light particles 

The spectrum of other light particles becomes also complicated in the presence of the R-symmetry breaking, since 
S a and M a scalars mix with each other via a cross term in the F -term potential \ms a S a — XSoM a \ 2 - To analyze the 
mass spectrum and interactions of those particles, we decompose M a and S a as 



S a = -^=(x s + iy s y a / f « 

M a = -^=(x m +iy m ). 
V2 



(26) 



Here, we have suppressed the index a, since the particles with different values of a decouple from each other in the 
following analysis. 

By using this expressions, we obtain a scalar potential, 



V = \X'A dyn M a 



^dyn — 2^ 



+ \\S a M a + \'A dyn S a \< 



\ ((A' 2 - A 2 )A 2 yn + ^(s + f R )^ xl + \ ((A' 2 + A 2 )A 2 yn 

1 1 AA' AA' 

+ ^ 2 ^lyn x l + 2 A ' 2A dyn2/s + -y|A dy „(s + f R ) x m x s + -^=A dyn (s + f R ) y m y s 

A 2 

+ ^ A dyn + 7^f( x m + Vm) ■ 



(27) 



Notice that the R-axion does not show up in the scalar interactions in this basis, and it only appears in the derivative 
couplings. From this potential, we find that the pseudo-Nambu-Goldstone mode resides not in (y m , y s ) but in (x m ,x s ). 
In the following, we concentrate on the real parts (x mi x s ). 
The mass-squared matrix of (x m , x s ) is given by, 



AA'Adyn (So) 

and hence, the masses of the eigenmodes (<fi, H) are; 



2 _ ( (A' 2 - A 2 )A 2 yn + A 2 (S f AA'A dy n (S ) 



m 2 H 



trM 2 
detM 2 



i (tr M 2 - v/(trAf 2 ) 2 - 4det 



M 2 ) = 



det M 2 



- (tr M 2 + y/(tr M 2 ) 



(2A' 2 - A 2 )A 2 yn + A 2 (S ) 
A' 2 (A' 2 -A 2 )A dy „ 



2 - 4 det M 2 
j ,< \ 2 



(28) 



(29) 

(30) 

(31) 
(32) 



The mixing angle is given by, 



Xm 



cost 
sin 6 



- sin 6 H, 
cos 9 H, 



(33) 



with 



tan 8 



AA'Adyn (So) 

A iv dyn " l 4 
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AA'A dyn (Sp) 
m\-m\ 



(34) 



As a result, we find that the lighter scalar cj> denotes the Nambu-Goldstone mode in the limit of A = A', and hence, 
we consider <j> as the dark matter. 

The R-axion interactions only appear in the kinetic terms. In the basis we have defined, the R-axion interaction 
comes from the kinetic terms of So and S a , 

C = l(da) 2 (l + -f) + -L(da) 2 (xl+yZ) + ^a(x s d»y s ~y s d»x s ). (35) 
* V JR/ 2 Jr JR 

Altogether, in the Nambu-Goldstone dark matter scenario (i.e., A' — A -C 1), light particles sector consists of the 
dark matter and the flaton in a TeV range, and the gravitino and the R-axion with much smaller masses, while the 
other components in S a and M a have masses of the order of the SUSY breaking scale. The most relevant terms for 
the dark matter annihilation is, then, given by, 

A 2 m 2 1 / \ 2 

C in t^-f R 5 2 s0 2 + -(da) 2 1 + -H . ( 36 ) 



fit 



where the first term comes from the scalar potential in Eq. (f2"T)) . while the second term comes from Eq. 



3. Flaton decay 

In order to discuss the dark matter annihilation via the s-channel exchange of the flaton, it is important to know 
the decay properties of the flaton. In particular, the decay rate into a dark matter pair is important even if the pole 
is unphysical, i.e., Im^ > m s , since the decay rate must be defined not on the exact pole, but on the center of mass 
energy of the dark matter collision, -Ecm- 

First, we consider the decay mode into a pair of the R-axions. The relevant interactions of the decay comes form 
the first term in Eq. ([36]) , and the decay rate into a pair of the R-axion is given by, 

1 TTl? 

62tt j R 

where we have neglected the mass of the R-axion and taken the final state velocity to be (3f = 1. For example, the 
decay rate is very small, i.e., T/m^z ICP 4 for /r > 30TeV and m 8 = 2TeV. As we will see this is favorable to realize 
a large effective boost factor via the Breit-Wigner enhancement. 

Next, we consider the flaton decay into a pair of the dark matter. The relevant interaction term is given in Eq. (|36p 
and the resultant decay rate is given by, 




r_^ = £E--ff , m„ (38) 



where /L denotes the size of the velocity of the dark matter. Notice that the value of Tg^^/pd, is well-defined even 
in the unphysical region, i.e., > m s . The value of T g^^/ (3<p is at most comparative to r s ^ aa , 

4 -3 



~ ( -— ) -f , (39) 



where we have used ~ m s /2 and vtlh S> m^>. Therefore, we find that the decay rate into a dark matter pair does 
not dominate over the one into an R-axion pair. 

Let us also consider the flaton decay into a pair of the gravitinos. The relevant interaction comes from the higher 
dimensional terms in the Kahler potential Eq. (|20[) , and the resultant interaction term is given by, 



Ant ~-yS# + / l .C. = ^ 5 # + / i .C, (40) 



= 

where we have used Eq. tfM$ , i.e., A 4 ~ F So /m s . Therefore, the decay width is suppressed by (to^/Asusy) 4 , an d 
hence, this mode is further suppressed compared with the mode into an R-axion pair. 
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Putting them all together, we obtain the flaton decay width at Ecm — m s , 

F s {Ecm) = r s ^ aa + Tg-,^ + • • • . (41) 

where Eqm > 2m0, and the dots refer to the modes into the MSSM particles (see appendix [C|) . In the following 
analysis, we approximate the above decay rate by, 

T S {E CM ) =s T s {m s ) ~ T„, (42) 

since all the other modes are subdominant at Ecm — m s . 

4- Dark matter annihilation via the s- channel flaton 

Now, let us consider the dark matter annihilation via the s-channel flaton exchange. The relevant interactions are 
again given in Eq. (|36|) . The amplitude of this process is given by 

M = X 2 2 ? M pm , (43) 

m H ~ m % e cm ~mj + im s r s (E C M) 

and the cross section by 



■Jrcl Pf ( ™% \ \ A E? 



aV Ie l 



CM 



32tt fo\m* H - ml ) (E CM -m^+ 



64tt \m H J ml (5 + + 7s 2 ' 



(44) 



where the T s and (3$ are defined at -Ecm > 2rri0. In the final expression, we have used the non-relativistic approxi- 
mation, 

E C m = + roji&t, (45) 

and introduced parameters S and 7 S by 

m 2 s =Aml{l-5), ls = T s /m s . (46) 

From this expression, we find that the annihilation cross section of the dark matter is substantially enhanced compared 
with the one given in Eq. (TIT)]) , for |<5|,7 S -C 1, which allows a sufficient annihilation cross section to reproduce the 
observed dark matter density. 

5. Dark matter density and Breit-Wigner enhancement 

Although we obtained the enhanced annihilation cross section of the dark matter, we should note that the thermal 
history of the dark matter density is drastically changed from the usual thermal relic density when the dark matter 
annihilates via the narrow resonance [24l [25| . and hence, the required annihilation cross section is different from 
the value given in Eq. Instead, in terms of the annihilation cross section at the zero temperature, the required 
annihilation cross section to obtain the correct abundance is given by , 

(™ re i)|T=o~lCr 9 GeV- 2 x f*. (47) 

x f 

Here Xf ~ 20 denotes the freeze-out parameter of the usual (non-resonant) thermal freeze-out history, while Xb is 
defined by, 

Xb (cUrel) \T=0 Jx f % 
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FIG. 1: Left) The <5 dependence of the required annihilation cross section at zero temperature from the observed dark matter 
density in the case of the unphysical pole. The red lines correspond to 7 S = 10 -3 , 7 S = 10 -4 and 7 S = 10 -5 from bottom to top. 
The green line shows the required annihilation cross section in the usual thermal history. The blue line shows the predicted 
annihilation cross section for A = 0.84, ran = 30TeV and m^> = 1 TeV. Right) The required annihilation cross section in the 
case of the physical pole. 



In the case of the unphysical pole, i.e., m s < 2m^, xj, is well approximated by min[<5 1 ,7 S ], and the above required 
annihilation cross section at the zero temperature is simply given by, 

(<">«*) |t=o ~ 1(T 9 GeV- 2 x — - . (49) 

Xf Max[d, 7 S J 

In Fig.[TJ we show the required annihilation cross section for given parameters as red lines. The figure shows that 
Eq. (|49p gives a good approximation. 

On the other hand, in the case of the physical pole, the estimation of Xb is much more complicated. In particular, 
the thermal average picks up the pole at v 2 cl — 4|<5| when the temperature is rather high, i.e., x~ x and hence, 

the annihilation cross section can be higher at the higher temperature than the one at the zero temperature. As a 
result, the required annihilation cross section at the zero temperature can be much lower than the one in the usual 
thermal relic history. In Fig.[TJ we also show the required cross section at the zero temperature in the case of the 
physical pole. The figure shows that the required cross section can be lower than the usual value. 

Now, let us compare these values with the dark matter annihilation cross section given in Eq. (|44| . For example, if 
we take, = 1 TeV, fn — 30 TeV, ran = 30 TeV and A = 1, the decay rate is very small 7 S ~ 10 -4 . In this case, 
the cross section at the zero temperature is 



(-01— ~-*x^ (l ^(^) 4 . (50, 

In Fig.[TJ we show the predicted annihilation cross section. From the figure, we find that the required annihilation 
cross section is achieved at 

5 ~ 10 -4 , (for unphysical pole), 

S ~ — 10 _1 , (for physical pole), (51) 

for the given parameter set (A = 1, = 1 TeV, and mu = fn = 30 TeV). Therefore, the Nambu-Goldstone dark 
matter is consistent with the observed dark matter density when it annihilates via the flaton resonance with the values 
of 8 given above. 

Interestingly, the Nambu-Goldstone dark matter predicts a non-trivial effective boost factor. The effective boost 
factor in the Breit-Wigner enhancement is defined by, 

BF = ^. (52) 

Thus, the effective boost factor for the above two solutions are given by, 

BF ~ 10 2 , (for unphysical pole), 

BF - 10~ 3 , (for physical pole), (53) 
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respectively for the parameters given above. 

Therefore, we find that the Nambu-Goldstone dark matter model predicts a non-trivial effective boost factor. 
Especially, the model with the unphysical flaton pole (m s < 2m is strongly favored in the light of the recent cosmic 
ray experiments. In this case, the parameter dependence of the boost factor is simply given by, 

Here, we have used Eqs. (|44]) and (|50|) . 



6. R-axion decay and anti-proton flux 

As we have discussed, the dark matter dominantly annihilates into an R-axion pair via the flaton resonance. 
Interestingly, since the R-axion has a mass in the range of tens to hundreds of MeV, it mainly decays into light lepton 
pairs (see Ref. |26| for detailed discussion on the R-axion properties). 

Therefore, the Nambu-Goldstone dark matter model provides a concrete example of the scenario developed in 
Ref. [l4j. There, the dark matter annihilate into a new light particle which subsequently decays into light leptons. In 
this way, we can obtain a hard positron spectrum without any additional anti-protons, so that we can explain the 
PAMELA results consistently. 

Here, we comment on the constraints on the decay constant and mass of the R-axion. For the R-axion in a mass 
range between two electrons and two muons, the stringent constraint comes from a beam-dump experiment [30(, which 
constrains the decay constant as 

fn > 10 4 5 GeV x ( Wa Y /2 , (55) 
JH ~ V 10 MeV/ ' V ' 

when we assume that the Higgs sector respects the R-symmetry and the R-charge of the so-called /z-term, H u Hd, is 
two [26| . Thus, our choice of the scales of the R-symmetry breaking and the SUSY breaking in the previous discussion 
are marginally consistent with the constraint. 4 On the other hand, for the R-axion heavier than two muons, the most 
stringent constraint comes from the rare decay of the T meson, Br(T — ► 7 + a) < 10~( 5 ~ 6 ) [3l|, which is given by 
fn > 10 3 GeV. Furthermore, since we are considering the R-axion with mass heavier than a few tens of MeV, it is 
free from the astrophysical constraints. 5 



IV. CONCLUSION 



In this paper, we have revisited the possibility of the dark matter in the SUSY breaking sector, in the light of the 
recent cosmic ray experiments. In our model, the dark matter is identified as a pseudo-Nambu-Goldstone mode in 
the SUSY breaking sector with a mass in a TeV range which makes it possible to interpret the observed bump in the 
e + + e~ flux at ATIC/PPB-BETS experiments. Interestingly, the observed dark matter density requires an existence 
of a narrow resonance through which the dark matter annihilates, which results in a large effective boost factor (in 
the case of the unphysical pole). In addition, the dominant final state of the annihilation process is a pair of the 
R-axions each of which decays into a pair of light leptons. Therefore, the Nambu-Goldstone dark matter model is 
quite favorable to explain the PAMELA anomaly. 

Several comments are in order. In the model of the Nambu-Goldstone dark matter, the SUSY breaking scale is 
around 30 TeV. Thus, the model is accompanied by the gravitino with a mass in a ten eV range. The gravitino with 
such a small mass is attractive, since it causes no problem in cosmology and astrophysics (27l. l28j . 

In the Nambu-Goldstone dark matter scenario, the dark matter mass is controlled by the degree of the explicit 
breaking of the approximate global SO (6) symmetry, i.e., the difference between A and A' (see Eq. (|ll|l ). One may 



4 For other choices of the R-charge of the /i-term, the constraint can be changed. For example, when the R-charge of the /i-term is zero, 
the R-axion does not mix with the neutral Higgs bosons in the SSM. In this case, the couplings between the R-axion and the SM fermion 
vanish at the tree-level, and hence, the above constraint on the decay constant is weakened. We may also consider the Higgs sector 
without the R-symmetry. In that case, the degree of the mixing between the R-axion and the Higgs bosons is also altered from the one 
discussed in Ref. 26] , which may weaken the above constraint. 

5 As discussed in Ref. :26) , the R-axion can be detected at the LHC experiment if the decay constant is in tens of TeV range which makes 
the R-axion mainly decay into a muon pair. 
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attribute the origin of the tuning between A and A' to a conformal dynamics at high energy scales. As discussed in 
Ref. [2!| , the conformal extensions of the vector-like SUSY breaking model possess an IR-fixed point where the global 
symmetry is enhanced. In such models, even if |A' — A| = 0(1) at a high energy scale, the couplings flow to the 
IR-fixed point in the course of the renormalization group evolution and end up with |A — A'| <C 1, at the scale of the 
SUSY breaking. Thus, if we assume that the SUSY breaking sector was in a conformal regime at higher energy scales 
than the SUSY breaking scale, we can explain the lightness of the dark matter compared with the scale of the SUSY 
breaking. 
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APPENDIX A: COLEMAN- WEINBERG POTENTIAL OF IYIT MODEL 



Here, we show the detailed analysis of the Coleman- Weinberg potential of the flaton Sq (see also Ref. [22].) The 
classical flat direction Sq is lifted by a one- loop correction via the interaction W — XSoM 2 . Using the notation 
(j = XSq, x = A' Adyn, y = A Adyn, the mass matrix for the fermions is 



and for the bosons 



M 2 = 



Mf = 


u 


o) 




x 2 + a 2 


— xa 


-y 2 





— X<7 


x 2 








-y 2 





x 2 + a 2 


— xa 








~xo 


x 2 



The eigenvalues of the fermion mass-squared matrix are 

m) = i (2a; 2 + cr 2 ± a^Ax 2 + cr 2 ) , 

while for bosons 

ml = ^ (2x 2 + a 2 -y 2 ± ^/\x 2 o 2 + cr 4 - 2y 2 a 2 + y 4 ) , 



2a; 2 + a 2 +y 2 ± ^J±x 2 o 2 + cr 4 + 2y 2 a 



2 -y 4 



Using this spectrum, we can compute the Coleman- Weinberg potential. 

5 



AVcw - 



64tt 2 



STr m 4 In m 2 



where a factor 5 comes from the number of M a . Expanding it up to second order in cr, we obtain 
STr m 4 mm 2 = -4a; 4 In a; + (x 2 - y 2 ) 2 ln(x 2 - y 2 ) + (x 2 + y 2 ) 2 \og(x 2 + y 2 ) 
+ 4 ((^ 2 + y 2 f log(a; 2 + y 2 ) - (a; 2 - y 2 ) 2 \og{x 2 - y 2 ) 

y 

-Ax 2 y 2 \og(x 2 ) - 2x 2 y 2 ) a 2 + 0(cr 4 ). 
Since y < x is needed to avoid tachyon, we take the small y limit as 

4w 4 

STr to 4 mm 2 = y 4 (3 + 41og(a;)) + ^ct 2 , 



(Al) 



(A2) 



(A3) 

(A4) 
(A5) 

(A6) 



(A7) 



(A8) 
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which is a good approximation even as y — > x. Within this approximation, the mass term for So from the Coleman- 
Weinberg potential is 



AVCW = ^3(\' A Cynf lXSl2 = 3^A^ A ^ |5 ° |2 - (A9) 

Therefore, we obtain the mass of the flat direction, 6 

mS ° ~ Vf (4^ Ady - (AU) 

Notice that the flat direction is also lifted by higher dimensional terms of So in the Kahler potential which is suppressed 
by the dynamical scale 47rAd yn - However, the flat direction mass is dominated by the one-loop contribution analyzed 
here, since the fields circulating in the loop is much lighter than 47rAd yn - 



APPENDIX B: R-BREAKING IN GAUGED IYIT MODEL 



In the main text, we have discussed how R-symmetry breaking in the hidden sector drastically changes the decay 
process of SSDM scenario. In this appendix, we study the R-symmetry breaking of the IYIT model with additional 
U(l) gauge symmetry. We embed 50(2) x 50(4) in the original 50(6) global symmetry of the IYIT model, where 
50(2) = 17(1) is gauged. The dark matter candidate M a lies in vector representation of 50(4) (a = 1, • • -4). 

The low-energy effective superpotential of the gauged IYIT model is given by 7 

W = xS+M- + xS-M + + yS a M a (Bl) 

with the constraint M+M- + ^M a M a — A 2 yn = 0. The subscript ± denotes the i7(l) charge of the chiral superfields. 
We parametrize the solution of the deformed moduli constraint as 

M+ = A dyn £>/V2Ad y „ t M _ = Adyn e -^/V2A dyn _ ( B2 ) 
In these variables, the leading order Kahler potential is canonically normalized: 

X=|5 + | 2 + |5-| 2 + |^| 2 +^(|M a | 2 + |5 a | 2 ) + -- - , (B3) 

a 

where non-canonical Kahler potential may be neglected when F$ ± = AA 2 yn <C A 2 yn . On the other hand, the 
superpotential can be written as 



W = mv (5+e-^/^ Ad - + 5_e^^ Ad -) ^/A 2 yn - -M a M a (B4) 

The tree level vacua have moduli space spanned by (f> = M a = S a = 0, 5+ = 5_ = a. At the tree level, massless 
degrees of freedom are one R-axion, one real modulus a = (Rc[5 + ] + Re[5_])/2 and one goldstino after gauging 
away the U(l) Nambu-Goldstone boson at generic points of the moduli space. The R-axion remains massless in the 
field theory limit, while the pseudo-modulus a will acquire a quantum potential, whose shape and resulting VEV 
determines whether the R-symmetry is broken. 

To see the R-symmetry breaking, we compute the Coleman- Weinberg potential V(a) = 
(1/647T 2 ) STrTO 4 (a) logm 2 (<r). In Fig. [21 we plot the Coleman- Weinberg potential for a given U(l) gauge cou- 
pling constant. The figure shows that the symmetry enhancement point a = becomes the local maximum and 
the Coleman- Weinberg potential develops a minimum at a ^ for a larger value of the additional gauge coupling 
constant. Thus, our U (1) gauged IYIT model serves as a perturbative model of R-breaking hidden sector with hidden 
dark matter. The R-breaking depends on the U(l) coupling constant. 



Corresponding Kahler potential corrections to reproduce AVcw ^ s given by 

5 A 4 



A,v - - wi^ (5 ° t5o)2 - (A10) 



7 When we set M a = S a = 0, our model is equivalent to k — > oo limit of the model with W = xS + M_ + xS-M + + kX(M + M- - ) 
which was studied in 1231. 
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FIG. 2: The Coleman- Weinberg potential for the U(l) gauged IYIT model. In the figure, we have assumed A = 0.5, A' = 1, 
and the gauge coupling constant as shown in the figure. 



APPENDIX C: F-FLATON DECAY INTO THE SSM PARTICLE 



In this appendix, we consider the decay modes of the fiaton into the SSM particles. Since we are assuming the 
model with gauge mediation, the fiaton couples to the SSM fields as the results of the mediation effects. For example, 
the effective coupling between the fiaton and the gauginos is given by a Yukawa interaction; 



-off 



1 m l 
2Jr~ 



1 + 



A 



s\ l X l + h.c. 



(CI) 



where mi denotes the gaugino mass and i runs the SSM gauge groups. Notice that the leading term in the above 
effective coupling is model independent as long as the messenger sector possesses the R-symmetry. On the other hand, 
the coupling between the fiaton and the sfcrmions depends on the messenger sector even if it is R-symmetric, and is 
given by, 



£eff = 



drrij 



ds 



X sff, 



(C2) 



s=0 



where the model dependent coefficient dm~/ds satisfies 



dm 2 ^ 



ds 



< 



s=0 



(C3) 



From these interactions, the fiaton decays into a pair of the SSM particles. For instance, the decay rate into a pair of 
the gluinos are given by 



r 



1 



4ir \m 



2 s 



(C4) 



Therefore, depending on the spectrum of the SSM and the dark matter (jtidm — m s /2), the branching ratio of the 
fiaton into the SSM particles can be suppressed. Notice that the branching ratio into the gravitino pair is highly 
suppressed [33} . 
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